#-*- coding: utf8 -*-
from tutor.scripts import *
from tutor.plugins.maple import *

# variáveis úteis
half = M('1/2')
t = M.t
_i, _j, _k = M('_i, _j, _k')
rr = M('`..`')
coeff = M.coeff

# parâmetros iniciais
a = M(oneof(1, 2, 3, 4, 6))
b = M(oneof(1, 2, 3, 4, 6))
c = M(oneof(1, 2, 3, 4, 6, 7))
D = M(oneof(1, 2, 3, 4, 6, 7))
Dy2 = D * y ** 2

# --- resultados --------------------------------------------------------------

# correto 
fluxo = a ** 2 * b ** 2 / 24
print('fluxo: %s' % fluxo)

# erros
fluxosup = fluxo + D * b ** 2 * c / 3
print('fluxosup: %s' % fluxosup)

stokes = D * b * c
print('stokes: %s' % stokes)

gauss = a * b * c / 120 * (3 * a ** 2 + 2 * b ** 2)
print('gauss: %s' % gauss)
